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Abstract 

A quantitative discussion of nonperturbative effects for the high tem- 
perature electroweak phase transition is presented. We propose a method 
for the computation of the temperature dependent effective scalar poten- 
tial that takes into account the running of the effective gauge coupling. 
Compared to perturbation theory we find a moderate decrease of the criti- 
cal temperature and an important change in the strength of the first order 
transition. We conclude that perturbation theory gives a misleading pic- 
ture of the dynamics of the transition. 
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1 Introduction 



It has recently been suggested (|T[ - ||) that the electroweak high temperature 
phase transition may be governed by a strongly interacting gauge theory. The 
associated nonperturbative phenomena are believed to play only a modest role 
if the mass of the Higgs scalar is very small. For realistic mass values consistent 
with the experimental lower bounds |J, however, the nonperturbative effects 
could dominate the quantitative and perhaps even qualitative behaviour of the 
high temperature effective potential for the Higgs scalar. If so, this will have 
important consequences for speculations that the observed baryon asymmetry 
might have been created during the electroweak phase transition in the early uni- 
verse. Baryon generation within the standard model not only requires substantial 
CP violation but also a sufficiently strong first order phase transition. This is 
necessary in order to guarantee a period in the evolution of the early universe 
where baryon number violating processes are out of thermal equilibrium. The 
order of the electroweak phase transition as well as many important details of its 
cosmological dynamics are encoded in the shape of the high temperature effective 
potential. This has recently been studied using different versions of resummed 
perturbation theory 0. Sizable deviations from the perturbative results change 
the rates of baryon number violating processes or bubble formation by many or- 
ders of magnitude and can therefore completely alter our picture of the phase 
transition. 

The importance of nonperturbative electroweak physics was first observed 
quantitatively §| in a simplified theory where all particles except for the Higgs 
scalar were neglected. In contrast to the perturbatively suggested || weak first 
order behaviour, the phase transition of this model was found to be of the second 
order, with critical exponents corresponding to the universality class of the three 
dimensional model. Near the transition the dimensionless couplings are strong 
(as dictated by the values of corresponding infrared fixpoints |TIJ ) and pertur- 
bation theory becomes inadequate. More generally, at nonvanishing temperature 
T the physics of the modes with momenta q 2 < (2ttT) 2 can always be described 
by an effective three dimensional theory fLT |. The crucial question is to what 
extent the fluctuations of these low momentum modes dominate the behaviour 
of the effective potential. In the pure scalar theory an investigation of the scale 
dependence of the effective average potential revealed that the running of the cou- 
plings becomes effectively three dimensional (||, for scales below 2nT, but 



the difference between the three- and the four- dimensional running has quantita- 
tively important effects only in a narrow temperature range around the critical 
temperature where the dimensionless quartic scalar coupling becomes strong. 

In the electroweak standard model one expects even stronger nonperturbative 
effects. In the effective three dimensional theory at scales below 2irT the running 
of the gauge coupling is not logarithmic as in four dimensions, but gouverned 
by a power law M. In consequence, the three dimensional "confinement scale" 
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was estimated between T/10 and T/5 ffl]. This scale is characteristic for the 
nonperturbative effects generated by the strong gauge coupling and determines, 
for example, the magnitude of W-boson condensates in complete analogy to QCD. 
In this work we give a more quantitative estimate of the modifications of the 
high temperature effective potential due to the three dimensional running of 
couplings. We will perform a "renormalization group improvement" for the scalar 
potential U(<p) in the effective three dimensional theory which is in spirit very 
similar to the "renormalization group improved one loop potential" proposed 
for the zero temperature case by Coleman and Weinberg ||13|| . This approach 
accounts properly for the three dimensional running of the effective couplings 
and is expected to be quantitatively very reliable as long as the (^-dependent 
gauge boson masses are sufficiently large compared to the confinement scale. 
This gives, in turn, a lower bound on the value of the scalar field \ip\ for which 
the potential U(<p) is quantitatively under control. For \ip\ smaller than this 
lower bound nonperturbative effects are expected to give important modifications. 
We will give rough estimates of the size of these effects and demonstrate that 
nonperturbative physics dominates the phase transition for realistic values of the 
mass of the Higgs scalar. 

The three dimensional physics is related to the infrared behaviour of the 
theory. Perturbation theory is plagued by strong infrared divergences in the 
presence of massless particles as the gauge bosons for ip = 0. A suitable method 
to deal with these difficulties is the average action r& |14| where an effective 
infrared cutoff k controls the infrared behaviour. The dependence of 1^ on the 
scale k is governed by an exact nonperturbative flow equation fL5 |. For nonzero 



temperature this flow equation accounts properly for the change from an effective 
four dimensional running of couplings for k 2 > (2ttT) 2 to the three dimensional 
running; for k 2 < (2ttT) 2 §. We use here a simplified description by a pure three 
dimensional theory for k < kx, &t = 2irT. The "initial value" of Tk T includes 
the effects of all quantum fluctuations except those of the modes of the three 
dimensional theory with q 2 < k\. It is obtained by integrating out the b / 
Matsubara frequencies, the high momentum modes of the n = frequency and 
all modes of the zero components of the gauge fields A Q . In particular we are 
interested in the average potential Uk and we can infer the shape of U\. T from 
the work of reference ]16|] (see below). For k — > the average potential becomes 
exactly the high temperature effective potential ([j§|,|15||). The latter can therefore 
be computed by solving the flow equation for k between kx and 0. As a further 
simplification we leave out the fermions and the hypercharge gauge field. 
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2 Nonperturbative flow equation 



Approximate solutions of the exact flow equation need a truncation of the most 
general form of T k and we will work with the ansatz 

T k [cp, A„] = J d d x (u k {p) + Z^ k \D^\ 2 + ^Z^F^F^ . (1) 

Here p = p^p and group indices are omitted. The form of the average potential 
U k (p) is left arbitrary and has to be determined by solving the flow equation. For 
the Abelian Higgs model the evolution equation for the average potential was 
computed in the approximation ([]]) in reference [jTj]]. Inserting the appropriate 
SU(2) group factors we obtain in the Landau gauge (a = 0), with t = Ink 

+ In (g 2 + fc 2 + m 2 ) + 31n(g 2 + fc 2 + m 2 )) (2) 
where the mass terms read 

m 2 B = \z^ k fp ; m\ = (U' k (p) + 2pU' k \p)) /Z^ k ; m\ = U' k (p)/Z^ k . (3) 

We have used here a masslike infrared cutoff R k = Z k k 2 but more general cutoff 
functions R k (q) may be employed. The partial derivative on the right hand 
side of (0) is meant to act only on R k and we omit contributions arising from the 
wave function renormalization Z k in R k . Primes denote derivatives with respect 
to p. 

This flow equation can be usedQ in arbitrary dimensions d and constitutes 
a nonlinear partial differential equation for the dependence of U on the two 
variables k and p. In our case it holds for the three dimensional potential Us and 
a correspondingly normalized scalar field p^. They are related to the usual four 
dimensional quantities by U 3 = U4/T , p 3 = p±/T. The gauge coupling in (Q) 
stands for the three dimensional running renormalized gauge coupling g 2 {k). Its 
value at the scale kx is given by 

*M = AMT (1 - ff ) (4) 

where 

ml = \gl{kr)T 2 (5) 

1 The ultraviolet divergence on the right hand side of equation (||) is particular to the use of 
a masslike infrared cutoff. For d < 4 it concerns only an irrelevant constant in U and is absent 
for dU'/dt. 
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and the effects of integrating out the A$ mode have been included in lowest order 



lq| . The evolution equation for the running gauge coupling in the pure Yang- 
Mills theory has been computed in reference ]T]. We use here mainly the lowest 
order result 

where the deviation of r from one accounts for the small contributions of scalar 
fluctuations which remain to be computed^. Furthermore, we need the anomalous 
dimension of the scalar field. For our purpose it can be approximated by ||T7 



* = J ~nt = ( 7) 



3 Running quartic coupling 

A convenient quantity for an investigation of the effective potential is the p- 
dependent quartic coupling 

V(p) = U»M = ^T" (8) 

Knowing for k = the function A 3 (p) = \3,o(p), the high temperature effective 
potential U(p) = Uq(p) can be reconstructed by integration and translation to a 
four dimensional normalization. One of the two integration constants is irrelevant 
and the other (the mass term linear in p) can be found by adapting U(p) to 
the perturbative result for large p where the three dimensional running of the 
couplings is irrelevant. In fact, for large p such that m B (p) > k\ one may use 
the one loop potential [TH 



U 3 (p 3 ) = - /i 2 (T)p 3 + ^(A 3 + AA 3 )p3 2 -Y^(6m| + 3m| + m? + 3^ (g) 
with ml = 3A 3 p 3 - p 2 (T), ml = A 3 p 3 - p 2 (T), Z v = 1. Here 

AAa - 64^r I, 1 + ^ 8^ Ml) (10) 

with Mh and M w the masses of the Higgs scalar and the gauge boson respectively 
is chosen such that (8ir 2 T 2 / gl(k T )) = A 3 . Two loop effects and corrections 

2 For a suitable choice of wave function renormalization constants in the infrared cutoff for 
the gauge bosons the lowest order result becomes independent of the gauge parameter a and 
can therefore be used for the Landau gauge employed in this paper. 
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from integrating out Aq may be included [IB]. We use here the mass term p 2 {T) 
as given by 

/AT) _ i Ml 37T _7T 
T 2 ~2T 2 4 " 4 M 2 1 J 

where a w = g\{hr)/^it, the weak fine structure constant at the scale kx- Here c = 
1 — -^p+Ac includes effects from integrating out the A mode. We use Ac = but 
higher loop corrections or the effects from including the quark fluctuations can be 
accounted for by an appropriate nonvanishing Ac. For example including the top 

quark and using 6 quark flavours for the Debye mass yields Ac = 3A ^ P + 7/^ U ^ 
where the second term accounts for the modification of txib- The initial conditions 
for the flow equation described below will be expressed for fixed p 2 {T)/T 2 and 
M|/M 2 . In consequence, Ac ^ results in a simple rescaling of temperature 
and the scalar field according to 

T~ od = T~ +^rAc 

<-Pmod/<P = T mod /T. (12) 

Here T mo d is appropriate for Ac ^ whereas T represents the scaling used in 
this work for Ac = 0. The above simple rescaling property is very useful for a 
quantitative comparison with authors using different prescriptions for the effective 
three dimensional theory. 

The evolution equation for the /c-dependence of A 3 (p 3 ) can be infered from 
(H) by differentiating twice with respect to p 3 and reads 

dXsAPs) 3 / 2% tk gl k k? 6Z- 2 fc A 2 ;fc (p 3 )A: 2 2Z- 2 A 2 , fc (p 3 )/c 2 \ 

dt 32n{ {k 2 +m 2 B fl + ( p + m2) 3/2 + ( p +m 2 )3 /2 ) ^ ) 

where we have neglected terms oc U^\p 3 ) and U^ 4 \p 3 ). We propose in this paper 
an approximate solution of the flow equation for A 3 (p 3 ) for k = 0. It is based on 
the observation that the p 3 -dependent mass terms m 2 B , m 2 , and m 2 act in equation 
(HD as independent infrared cutoffs in just the same way as k 2 . A variation of 
m 2 for k 2 = is roughly equivalent to a variation of k 2 at m 2 = 0. We use 
this observation to translate the flow equation (p~3|) into a renormalization group 



equation for A 3 (p 3 ) at k = 0: In equation flTB| ) we replace by -j^ = niB-g^ and 



the factors k 2 ik 2 + m 2 ) 3 ^ 2 by rrT 1 . In order to see that our simple prescription 
indeed corresponds to an approximate solution of the flow equation ( f[3|) we omit 
for a moment the two last terms in (|i~3|) which arise from the scalar fluctuations. 
Writing 

I (™»sb^<*>) = ^'^s ( fc2 + m| ) " /2 (14) 
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we observe that the right hand side is suppressed both for k m B and k <C m B . 
Since the evolution of the quantity m B Q^Xs is dominated by a narrow interval 
k ~ m B , we may approximate on the right hand side of equation (14] ) g%(k) and 
Ztp(k) by g^(m B ) and Z^wb), and similar in the relation (|) between and 
p 3 . Equation ( |14D is then easily integrated between = and k = k T to give 



9 3 (9 

m BT{ A 3 , (p 3 ) = —Z 2 (m B )gl(m B )m B l + m B - A 3 , fcT (p 3 ) - 

om B 327T ^ " om B 

— —Zl(m B )gUm B ) m * , /0 . (15) 

The sum of the last two terms is small compared to the first termQ which coincides 
exactly with our prescription. This type of arguments can be generalized for the 
scalar fluctuations. In summary, we can now work with a new effective infrared 
cutoff k' = m B which is a function of p 3 (we omit the prime on k in the following) , 



1 



m. 



-Z v (k)g!(k)p 3 . (16) 
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This procedure transforms equation flT3| ) into a simple differential equation for 
A 3 (p 3 ) = A 3 (&(p 3 )). m terms of the renormalized coupling 



X R (k) = Z~\k)Uk) (17) 



it reads 



m^ k) = 3^k ^ {k) + + ^) ^ /2 (^3 - ^9l(k)X R (k). (18) 

For the terms oc — and oc ^ from equation (|15D we have approximated in (|3|) 
£/'(p 3 ) ~ p 3 A^(p 3 ) which amounts to neglecting the mass term. For negative 
X B (k) our approximation does not describe properly the effect of the scalar fluc- 
tuations. Since their contribution is small in this region we simply omit the terms 
oc X 3 B 2 once X R (k) becomes negative. 

The renormalization group equation ( |T8"D for the p 3 - dependence of U" is the 
central equation of this work. Except for the last term arising from the anomalous 
dimension it can be directly obtained by taking appropriate derivatives of the 
one loop formula (§), treating mass ratios such as m B jm\ as k- independent and 
replacing at the end the couplings g\ and A 3 by running couplings evaluated at 
the scale k. For Z v = 1 and g 2 , A 3 independent of k equation (|18D reproduces 
exactly the one loop result. Our renormalization group equation enables us to 

3 This statement is immediately apparent for tob <C fcr where both terms are seperately 
small. We expect it to hold also for tub ~ fcr, but a proof needs a detailed understanding of 
the transition region between the effective three dimensional and the effective four dimensional 
theory. 
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include the effects of running couplings and the anomalous dimension. Combining 
equations @, (§), (0), ([H]) and ([18]) we can compute the p 3 -dependence of the 
high temperature effective potential by a solution of the flow equation. The initial 
values for this solution are set by equations (H), (|TT|) and 



where the last value is extracted from the validity of the one loop potential (||) 
at fcr- 



4 Analytical approximations 

Before turning to the numerical solution of the flow equations we display some 
analytical results in appropriate approximations. 

First, from (^) one easily finds that the running gauge coupling is given by 

1 1 23r / 1 1\ 

+ 7T7- T--T • (20) 



gl(k) gi(hr) 24vr \k T k 
Thus, g%(k) diverges as k approaches koo with 



, 1 24vr > 



This reflects the well known fact that the SU(2) Higgs model is confining in 3 
dimensions. From (|7]) we can compute Z v (k) and find 



6 

" 23t 



Z M = ( ) » ( 22 ) 



T 



where we set Z^kr) = 1- Using (|22|) we can explicitly relate the renormalized to 
the unrenormalized couplings. Also, ps(k) is given by 

k 2 ( 23ra w ( 2txT\\ x -^ . , 

*w = w( 1 + ^r( 1 --r)) ■ (23) 

For small Ar(/c) we can neglect the term oc A^ 2 on the right hand side of (|i~8|) 
and solve this equation analytically. One easily recovers the 1-loop result if one 
sets rjip = and g\ = const.. We can also find a solution if we include r} v and the 
running of the gauge coupling. Defining 



S 



the flow equation reads 

dR 9/12 



dlagl 92r V 23r 
This is easily solved for R and one finds 



R. (25) 



A «Miif - ra^ 1 - (If i i- < 26 > 

The main changes as compared to the 1-loop calculations can be understood 
from the corresponding differential equations: The inclusion of 7] v lowers the 
scale at which U" changes sign, whereas the running of the gauge coupling acts 
the opposite way. Thus at a given k, corresponding to a given p 3 , the running 
of g\ makes the potential bend up less than the 1-loop calculation predicts. The 
inclusion of the neglected terms oc Xr 2 will have the same effect as soon as they 
are not small, i.e. for larger Higgs masses. 

In fig. 1 we have depicted the quartic coupling U"(p) (in a four dimensional 

1 /2 

normalization) as a function of $ = p 4 for Mj, = 35 GeV. The corrections due 
to the running gauge coupling and the anomalous dimension are small for $ > 
100 GeV. There are sizable modifications in the region 40 GeV < $ < 60 GeV 
and large corrections for $ < 40 GeV. Translated to the effective potential the 
running of g\ should strengthen the phase transition and lower the critical tem- 
perature. This is what one would naively expect, since the first order character 
of the transition is due to the gauge-boson loops. Thus, enhancing the coupling 
should give a transition more strongly first order. We should note, however, that 
the anomalous dimension has the opposite effect. To estimate the precise effects 
of these modifications on the shape of the potential we have to proceed to a 
numerical integration of equation 



5 Numerical solution of the flow equation 

In order to make the numerical solution as easy as possible, we have transformed 
the system of differential equations (§) and (|T8| ) into a third order differential 
equation for U 3 . Since equation ( ^3|) gives p 3 in terms of k and is not easily 
inverted analytically, we use k as the integration variable. This yields 

d 3 U 3 „ 2 2 „ m Bd 2 U 3 ( B 2 C\dU 3 . . 
3 - = A 2 Z 2 ^ + 3-—^ - 3— - - — i (27) 



dk 3 ^ a dk 2 



where 



32nk 2 
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As initial conditions we use 



P&TJ f)TJ 

^r(fcr) = A 2 {k T )\ 3 {k T ) + B(k T )^(k T ) 
~^{k T ) = A(k T )^(k T ) ; U 3 (k T ) = U 3 (k T ) ; Z v {k T ) = 1. (29) 
Us represents the 1-loop potential (§) such that 

= V(r) + (A3+AA 3 )^-^T^I(^) + (V6+^)^ T ) .(30) 



T 



A similar set of equations relates the derivatives of the potential with respect 
to p 3 to the solution of (p7|) . We integrate equation (p7|) with the help of the 



Runge-Kutta method and use for all numerical work g 4 = 2/3, M TO = 80.6 GeV, 
and t — 1. 

In fig. 2 we show the effective potential as obtained with our method for dif- 
ferent temperatures and for masses of the Higgs scalar of 35 and 80 GeV. In 
all the plots we use four dimensional quantities and plot AU = U($) — U(0). 
For M/j = 35 GeV (fig. 2a) the first two curves correspond to the critical tem- 
perature as given by the one loop approximation (T = 97.50 GeV) and the crit- 
ical temperature obtained from our renormalization group improved approach 
(T = 95.85 GeV). Even though the critical temperature is not changed by much, 
the shape of the potential at a given temperature varies considerably. At the one 
loop critical temperature there remains not even a local minimum once renormal- 
ization group effects are included! We also find important quantitative differences 
for the size of the barrier at the critical temperature in both approaches. This 
effect becomes larger for higher values of the Higgs mass. For Mh = 80 GeV 
(fig. 2b) the critical temperatures (185.23 GeV for one loop, 184.25 GeV for RG- 
improvement) remain similar, but the shape of the potential differs strongly 
(even though this is difficult to see on the scale we use in fig. 2b) . We have 
also indicated by a square the value of <3> where the dimensionless gauge coupling 
g 2 (k) = g 2 {k)/k becomes as large as 27r, and similarly, by the end of the solid line, 
where it reaches Atc. These points may be considered as an estimate for the limit 
of validity of our approximations. To the left of this region we have to deal with 
a strongly interacting gauge theory. Straightforward use of our approach in this 
region would yield the potential as indicated by the diamonds. We observe that 
our renormalization group improved effective potential can be formally extended 
to $ = 0: The confinement scale k^ corresponds to p = since by equation (|i~6|) 
p 3 oc g 3 2 (k)k 2 and g 3 {koa) = oo. This explains the different qualitative behaviour 
as compared to |J. 

In order to demonstrate that g 2 ~ 2tt — Arc indeed corresponds to the onset 
of the strong coupling regime we notice that g 2 (k) = 2tt is reached at a scale 
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k ~ 1.5/cqo, which is already very close to the confinement scale (cf. eq. (p0|) ) . 
The anomalous dimension r] v reaches the value —1 for g 2 (k) ~ Arc (cf. eq. (^)). 
Instead of the lowest order estimate for the running of g 2 (|6|) we may also use the 
improved nonperturbative estimate of reference ffl, 



dgi(k) 23 



dt 24nk \1 - §^gi(k)k- 



f 3 (k). (31) 



Here the /3-function diverges^ for g 2 {k) = ||7r, which is a value close to 2n. For 



the /3-function (RTI) the confinement scale turns out higher and corresponds now 
to a nonzero value of $. We compare the results of the /3-functions (||) and (|3T|) 
for the effective potential in fig. 3, where we also display the one loop potential 
at the same temperature. The diamonds correspond to fl3T|) and this curve ends 
once (3 g 2 diverges. For ease of comparison we have shifted all three curves to 
conicide at $ = 100 GeV. The temperature chosen is the critical temperature for 
Mh = 80 GeV as obtained by the renormalization group improved potential with 
the lowest order /3-function ([]). We emphasize that strong interaction phenomena 
cover here the whole region between the two minima. Any perturbative estimate 
of details of the potential as needed for calculations of bubble nucleation and for 
a treatment of the cosmological dynamics of the phase transition seems highly 
unreliable for Higgs masses of about 80 GeV or higher. 

At this point we should mention that the renormalization group improvement 
has already extended considerably the range of $ where a reliable computation is 
available: In the standard perturbative loop calculations the 92-dependence of the 
effective gauge coupling is usually not taken into account. Since the gauge-boson 
fluctuations are crucial for creating the barrier responsible for the first order 
transition, a change of g 2 by a factor of two as compared to the perturbative 
value (which in our language corresponds to g 2 (kx)) seems to be at the limit of 
what is tolerable for a perturbative estimate to be quantitatively correct. From 
( |2"0"D we find that g 2 (k) j 'g 2 {k,T) = 2 corresponds to a value of $ = 0.55T (we have 
indicated this value by a cross on the curves in fig. 2). In contrast, the value 
where g 2 (k) = 2tt (squares in fig. 2) corresponds to <3> nj) = 0.28T which is smaller 
by a factor of two. Inspection of figures 1 and 2a suggests that even the range 
$ > 0.55T for the validity of the loop expansion is overestimated and a more 
conservative estimate amounts to $ ~ T. 



6 Strong electroweak interactions 

We finally turn to the region of strong interactions to the left of the solid line 
in fig. 2. We first observe that our computation of the potential has effectively 

4 This is, of course, not a physical effect; see the discussion in Jf] 
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included only the quantum fluctuations with q 2 > k 2 ^ since even for $ = the 
infrared cutoff is given by k^. In consequence, the potential we have computed 
and shown in figures 2 and 3 corresponds to the average potential Uk x rather 
than to the effective potential Uq. For $ sufficiently large, i.e. for $ > $ np to the 
right of the squares in fig. 2, the difference between Uq and Uk x should be small 
since the contribution of the low momentum modes is suppressed by an effective 
infrared cutoff rug ^> k^ . This is not true anymore in the region of small $ where 
the low momentum fluctuations are expected to produce large nonperturbative 
effects. In particular, we notice that the ansatz ([]]) becomes insufficient: At the 
scale fcoo the gauge boson kinetic term oc Zpik,^ vanishes Zp(k) and becomes 
negative |J for k < k^. This is a clear indication that for k = the minimum 
of the Euclidian effective action does not occur for F llv F tlv = but rather for 



a nonzero value of F 2 |19| . One expects W-boson condensates in close analogy 



to the gluon condensates in QCD. The best way to visualize these effects in the 
context of the effective average action is perhaps the introduction of a composite 
scalar field for the operator ^F^F^. This can be done at a scale in the 
vincinity of (somewhat above) k^ according to the general formalism proposed in 
[pO|1 . The average scalar potential is then generalized to a potential depending on 
two scalar degrees of freedom, O). For k > k^ the minimum of Uk(<&, O) 

occursf] at < >= whereas for k — > the minimum value < > will be a 
nonvanishing function of $ with <0> ($) > for $ smaller than some critical 
$ cr . In view of fig. 3 we may roughly associate the critical $ cr for the onset of 
condensation phenomena with $ np , denoted by the squares in fig 2. The W-boson 
condensate will lower the value of the effective potential U(&) = U(Q, <0> ($)) 
in the region of small <3> as compared to the computed U^^jO) [0. We will 
present in the following some quantitative but rough estimates how this affects 
our picture of the phase transition. 
Let us try to estimate the difference 

<y 1 = y fcoo (o,o)-E/ (o,<e>(o)) (32) 

which measures how far the W-boson condensation lowers the potential at the ori- 
gin. Since for $ = the only relevant scale for the condensate is the confinement 
scale fcoo, we can infer from simple dimensional analysis 

AU 3 = Kkl (33) 

with K a constant of order one. Since the running of the gauge coupling in three 
dimensions obeys a power law (in contrast to the four dimensional logarithms), 
the scale k^ is determined relatively precisely once the complete beta- function is 
known. In the relevant region the gauge coupling is large and there is no apparent 

5 Strictly speaking <0> may take a nonzero value depending on its precise definition. For 
the purpose of this qualitative discussion we shall speak of <&>= if the effect of <6> is 
small. 
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small dimensionless parameter in the problem and we will therefore use K = 1. 
With fcoo = 0.1 3T we infer for Si 

61 = 2.2 x 10~ 3 ifT 4 . (34) 

We have depicted Si in figure 2 and also show the potential for the "nonpertur- 
bative critical temperature" corresponding to T = 94.70 GeV(173.74 GeV) for 
M h = 35 GeV(80 GeV). For M h = 35 GeV the change in the critical tempera- 
ture is not very large, but the size of the barrier as obtained by extrapolating 
the potential from $ = $ rap to the value AU = —S\ for $ = is considerably 
enhanced by the W-boson condensation. For Mh = 80 GeV this effect is even 
more dramatic, and also the critical temperature is lowered considerably. A com- 
bination of our renormalization group improved potential for $ > $ np and other 
nonperturbative estimates of the critical temperatures using methods sensible to 
condensation phenomena - for example lattice calculations - can be used to de- 
termine the proportionality factor K. For this purpose we plot in figure 4 the 
critical temperature as a function of the Higgs-mass for various values of K. Here 
we have also included the (properly rescaled) results of recent lattice studies for 



M h = 35 GeV [^TJ and 80 GeV J^l respectively. These values seem to support 
the general picture, yielding K of the order of 10. Conclusive evidence surely 
needs more data however. 

It is also interesting to obtain an estimate for a lower bound on the critical 
temperature. For this purpose we use 

S 2 = U koo ($ np , 0) - U (0, <6> (0)) = K'kl p T ~ 8 x 10- 3 iTT 4 . (35) 

Since, as argued above, k np denotes the scale for the onset of strong interaction 
phenomena it is difficult to conceive that the proportionality constant can exceed 
one by much. We display S2 for K' — 1 together with the potential for the 
corresponding critical temperature in fig. 2. The "lower bounds" on the critical 
temperature are 

T cr > 92.54 GeV for M h = 35 GeV 

T cr > 159.14 GeV for M h = 80 GeV (36) 

and somewhat below these values for K' > 1 . It is apparent from fig. 2 that 
the lower bounds on T cr also correspond to upper bounds on the surface under 
the barrier which determines the strength of the first order phase transition. We 
observe that the use of the improved /3-function ([H|) enhances fcoo by a factor of 
two and leads to the guess K = 8. This yields T cr in the vincinity of the "bound" 
(|36D, and seems to be favoured by the lattice results available. 

Our estimates can be extended to larger values of the Higgs mass which are up 
to now difficult to access by alternative methods. In fig. 5 we plot the effective po- 
tential for Mh = 140 GeV and 200 GeV for three values of the temperature: The 
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upper curve corresponds to the critical temperature as determined from in 
the absence of condensation phenomena, the middle curve shows our estimate us- 
ing K = 1 for the critical temperature including W-boson condensation, and the 
lower curve corresponds to the "lower bound" given above. For Mh = 200 GeV 
we notice that $ np (again indicated by a square) almost approaches the mini- 
mum of U for T cr = 318.9 GeV. This may have important consequences for the 
dynamics of the phase transition: It is conceivable that there is almost no barrier 
between the minima at $ = and $ ^ 0, and that the barrier is even completely 
absent for Mh ~ 200 GeV. This could be interpreted as a change from a first 
order transition to an analytical crossover for very large Higgs masses [0. Even 
though there is at presence no evidence for an analytical crossover, this possibility 
cannot be excluded for Mh ~ 200 GeV! 



7 Conclusions 

In conclusion we have presented here a renormalization group improved estimate 
of the high temperature effective potential which determines the dynamics of the 
electroweak phase transition in the early universe. It is based on an approximative 
solution of an exact nonperturbative flow equation and includes properly several 
effects not accounted for in the presently available results from the loop expan- 
sion. This concerns, in particular, the running of the effective three dimensional 
gauge coupling and the anomalous dimension of the scalar field. We believe our 
estimates to be quantitatively accurate for large enough values of the scalar field, 
$ > $„ p . In this region of field space, the differences as compared to the loop 
expansion are already considerable even for a mass of the Higgs scalar as low as 
35 GeV, being further enhanced for large scalar masses. We argue that nonper- 
turbative effects not included in the present calculation further lower the critical 
temperature. They also enhance the barrier characteristic for the strength of the 
first order phase transition for low and moderate scalar masses {M h ~ 100 GeV). 
This may lead to a sufficient strength of the first order transition to be compatible 
with electroweak baryogenesis for realistic masses of the Higgs scalar. For large 
scalar masses (Mh > 200 GeV) it is not excluded that the barrier disappears and 
the phase transition turns into an analytical crossover. We have also estimated 
lower bounds on the critical temperature including nonperturbative effects. The 
main uncertainty in the present calculation concerns the size of nonperturbative 
condensates (W-boson condensation), which determine the behaviour of the po- 
tential for small values of the field. In view of the extreme sensibility of tunneling 
rates to the height of the barrier and the crucial importance of nonperturbative 
phenomena for an estimate of this barrier, we believe that a quantitative under- 
standing of the dynamics of the electroweak phase transition has to wait for a 
reliable treatment of condensation phenomena by nonperturbative methods. 
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Figure Caption 



Fig. 1 Different approximations to the effective quartic scalar coupling for Mh = 
35 GeV Higgs at the RG-improved critical temperature. 

Fig. 2 Renormalization group improved effective potential for Mh = 35 and 80 
GeV at different temperatures. 

Fig. 3 Comparison of the effective potential using different /5-functions for the 
gauge coupling. 

Fig. 4 Critical temperature including condensation effects as a function of the 
Higgs-mass for different values of K. 

Fig. 5 Renormalization group improved effective potential for M h = 140 and 200 
GeV at different temperatures. 
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